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Summary 


A  theoretical  linear  stability  analysis  is  used  to  consider  the  effect  of  a  porous  wall  on  the  first 
mode  of  a  hypersonic  boundary  layer  on  a  sharp  slender  cone.  The  effect  of  curvature  and  of 
the  attached  shock  are  included  for  axisymmetric  and  non-axisymmetric  disturbances.  The 
flow  in  the  hypersonic  boundary  layer  is  coupled  to  the  flow  in  the  porous  layer.  The  linear 
results  for  neutral  stability  and  spatial  stability  are  presented  for  physical  parameters  and 
porous  wall  models,  chosen  to  correspond  to  relevant  experiments.  The  effects  of  varying 
the  porous  wall  parameters  are  investigated.  A  weakly  nonlinear  stability  analysis  is  carried 
out  allowing  an  equation  for  the  amplitude  of  disturbances  to  be  derived.  The  coefficients  of 
the  terms  in  the  amplitude  equation  are  evaluated  for  axisymmetric  and  non-axisymmetric 
disturbances.  The  stabilising  or  destabilising  effect  of  nonlinearity  is  found  to  depend  on  the 
cone  radius.  The  presence  of  porous  walls  significantly  influences  the  effect  of  nonlinearity. 
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1  Introduction 


Transition  to  turbulence  in  hypersonic  flows  is  associated  with  amplification  of  the  first 
and/or  second  Mack  modes.  The  first  Mack  mode  is  the  high  speed  counterpart  of  Tollmien- 
Schlichting  waves,  so  a  viscous  instability,  with  modes  located  close  to  the  boundary.  The 
second  Mack  mode  is  an  inviscid  instability.  The  second  Mack  mode  is  believed  to  be 
responsible  for  transition  to  turbulence  on  hypersonic  slender  bodies.  Recent  experiments 
by  Fedorov  et  al.  [1-3]  have  shown  that  a  porous  coating  greatly  stabilizes  the  second  mode 
of  the  hypersonic  boundary  layer  on  sharp  slender  cones.  The  effect  of  the  porous  coating 
is  to  reduce  the  growth  rates  of  the  second  mode  to  a  level  where  they  are  comparable 
with  those  of  the  first  mode  (occurring  at  lower  frequencies).  In  addition,  the  first  mode  is 
observed  to  be  slightly  destabilized  by  the  presence  of  the  porous  coating.  Thus,  the  first 
mode  may  now  be  more  significant  in  the  transition  process. 

Previous  work  [4]  on  the  linear  instability  of  the  first  modes  (associated  with  Tollmien- 
Schlichting  waves)  in  the  hypersonic  flow  over  a  cone  includes  the  effect  of  curvature  and  also 
the  attached  shock.  The  effect  of  curvature  is  significant  with  disturbances  only  existing  over 
a  finite  range  of  wavenumbers  for  a  fixed  radius.  Beyond  a  critical  radius,  dependent  on  the 
azimuthal  wavenumber,  all  disturbances  are  damped.  In  [4]  it  was  demonstrated  that  the 
effect  of  the  shock,  which  allows  incoming  and  outgoing  waves,  gives  rise  to  multiple  modes. 
It  was  demonstrated  that  the  effect  of  the  shock  should  not  be  neglected  otherwise  the  correct 
effect  of  curvature  will  not  be  realized.  The  modes  which  exist  in  the  absence  of  the  shock  are 
now  totally  destroyed.  The  effect  of  nonlinearity  and  curvature  in  hypersonic  boundary  layer 
flow  has  been  considered  in  [5].  The  current  study  investigates  the  effects  of  a  porous  layer 
on  the  first  mode  (viscous)  for  axisymmetric  and  non-axisymmetric  disturbances.  This  is  to 
confirm  the  experimental  observations  with  a  theory  which  includes  the  effect  of  curvature. 
The  effect  of  curvature  is  of  practical  importance  and  has  not  been  previously  investigated 
theoretically  for  a  porous  coating.  Also,  if  the  second  mode  (and  higher  modes)  are  greatly 
stabilized  then  the  first  mode  may  well  lead  to  transition  to  turbulence,  particularly  if  the 
disturbances  are  triggered  by  wall  roughness,  since  this  mode  is  located  close  to  the  wall  and 
governed  by  viscous  effects. 

We  consider  the  effect  of  a  porous  wall  on  the  linear  instability  and  weakly  nonlinear 
stability  of  hypersonic  flow  over  a  sharp  slender  cone.  In  this  theoretical  and  asymptotic 
investigation  for  large  Mach  number  and  large  Reynolds  number  the  scales  used  will  be 
appropriate  to  the  first  mode  instability  which  is  governed  by  a  triple-deck  structure.  The 
effects  of  curvature  and  the  attached  shock  will  be  taken  into  account.  The  effect  of  the 
porous  wall  will  change  the  boundary  condition  on  the  normal  velocity  at  the  interface. 

Our  previous  results  have  considered  the  effects  of  three  different  porous  wall  models 
used  by  other  investigations;  the  porous  wall  model  of  Fedorov  et  al.  [1]  for  a  regular 
microstructure  and  a  random  microstructure  [2]  and  a  mesh  microstructure  used  in  the 
experiments  of  Lukashevich  et  al.  [6].  The  latter  authors  considered  hypersonic  flow  over 
layers  of  stainless  steel  wire  mesh.  It  was  demonstrated  experimentally  and  theoretically 
that  a  thinner  layer  of  this  structure  is  required  to  give  optimal  stabilisation  of  the  second 
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Mack  mode.  Thus,  it  may  be  possible  that  this  type  of  regular  microstructure  could  be 
incorporated  into  thermal  protection  systems. 

These  porous  wall  cases  are  characterized  by  an  admittance  Ay  which  is  a  function  of 
the  disturbance  frequency  and  depends  on  the  physical  properties  of  the  flow  and  the  porous 
layer.  The  focus  of  this  report  is  the  effect  of  these  types  of  porous  layers  on  the  nonlinear 
stability  of  the  first  Mack  mode. 

In  the  first  report  the  linear  results  were  presented  for  axisymmetric  disturbances.  It 
was  demonstrated  that  the  effect  of  the  porous  wall  was  to  reduce  the  neutral  values  slightly 
and  gave  rise  to  significantly  larger  growth  rates  than  those  for  a  solid  wall  for  some  pa¬ 
rameter  values.  The  second  report  considered  linear  stability  for  non-axisymmetric  distur¬ 
bances.  It  was  shown  that  the  destabilising  effect  of  the  porous  wall  was  more  significant  for 
non-axisymmetric  disturbances.  The  third  report  considered  the  effects  of  the  porous  wall 
properties  on  the  maximum  spatial  growth  rates  obtained  from  the  linear  stability  analysis. 
In  particular,  the  effects  of  increasing  pore  radius,  porosity  and  pore  depth  were  shown  to 
increase  the  spatial  growth  rates,  although  they  did  level  off.  The  effects  of  gas  rarefaction 
were  also  investigated.  In  addition,  a  nonlinear  analysis  was  presented  which  lead  to  the 
derivation  of  an  amplitude  equation.  The  fourth  report  considered  the  linear  stability  of  the 
hypersonic  flow  over  a  random  microstructure  using  the  theoretical  model  of  Fedorov  et  al. 
[2]  appropriate  to  a  felt  metal  surface.  It  was  found  that  the  random  microstructure  has  a 
larger  destabilising  effect  than  the  regular  microstructure,  particularly  for  non-axisymmetric 
disturbances.  Much  larger  growth  rates  were  obtained  for  the  random  microstructure  than 
for  the  regular  microstructure.  Thus,  in  a  practical  setting  it  is  expected  that  a  random 
microstructure  may  cause  transition  to  turbulence  if  first  modes  are  excited.  These  results 
agree  with  the  numerical  simulations  of  Wang  and  Zhong  [7]  who  compared  the  effect  of 
a  felt-metal  porous  coating  with  that  of  a  regular  structure  porous  coating  on  Mack’s  first 
and  second  modes.  The  fifth  report  considered  the  flow  over  layers  of  fine  mesh  of  square 
cross-section,  motivated  by  the  recent  experiments  of  Lukashevich  et  al.  [6],  The  sixth 
report  presented  some  results  from  the  weakly  nonlinear  analysis  for  non-axisymmetric  dis¬ 
turbances.  These  results  showed  that  the  nonlinear  effects  were  destabilising. 

This  final  report  summarises  all  our  findings.  For  simplicity,  we  consider  flow  relating  to 
the  previous  experimental  studies  [1,2]  and  [6].  Thus,  we  concentrate  on  presenting  linear 
and  nonlinear  results  relating  to  physical  situations. 

In  §2  the  problem  formulation  is  given  for  flow  over  the  three  types  of  microstructure 
previously  investigated.  In  §3  the  linear  stability  problem  is  summarised  for  axisymmetric 
and  non-axisymmetric  disturbances,  resulting  in  two  dispersion  relations.  §4  presents  the 
results  corresponding  to  experiments  for  neutrally  stable  disturbances,  while  §5  considers 
the  effect  of  the  porous  walls  on  the  spatial  growth  rates.  In  §6  the  effect  of  varying  the 
porosity,  pore  radius  and  pore  depth  are  illustrated  by  considering  the  maximum  spatial 
growth  rates.  Comparisons  are  also  made  for  a  porous  surface  optimised  for  second  Mack 
mode  stabilisation.  In  §7  the  nonlinear  analysis  is  outlined  and  the  results  are  presented  for 
axisymmetric  and  non-axisymmetric  disturbances  in  §8.  In  §9  we  draw  some  conclusions  of 
our  study. 
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Shock 


Figure  1:  The  geometry  of  the  cone  and  shock.  The  cone  is  taken  to  be  of  semi-angle  9C  with 
the  attached  conical  shock  making  an  angle  9S  with  the  surface  of  the  cone. 


2  Methods,  Assumptions,  and  Procedures 


The  flow  of  a  compressible,  viscous  fluid  over  a  sharp  cone  with  a  porous  boundary,  of  semi¬ 
angle  9C  is  considered  at  hypersonic  speed  Uq  aligned  with  its  axis.  We  consider  an  attached 
conical  shock  which  makes  an  angle  9S  with  the  cone;  a  situation  which  is  illustrated  in  figure 
1.  Spherical  polars  ( x ,  9,  <fi)  is  the  natural  coordinate  system  in  which  to  describe  the  basic 
flow,  and  here  (j)  denotes  the  azimuthal  angle.  Furthermore,  the  radial  distance  x  has  been 
non-dimensionalised  with  respect  to  L*,  the  distance  from  the  tip  of  the  cone  to  the  location 
under  consideration. 


The  approximate  basic  flow  used  is  described  in  [4]  so  the  complete  details  are  omitted. 
The  important  features  are  summarised  below.  Away  from  the  surface  of  the  cone  the  flow 
satisfies  the  (inviscid)  Euler  equations.  The  velocities  are  non-dimensionalised  with  respect 
to  [/_,  where  f7_  is  the  magnitude  of  the  fluid  velocity  just  behind  the  shock.  Additionally, 
the  time,  pressure  and  density  are  non-dimensionalised  with  respect  to  L*/U-,  p-U'i  and  p_ 
respectively,  where  p_  is  the  density  just  behind  the  shock.  Finally,  the  basic  temperature 
is  non-dimensionalised  by  T_,  the  temperature  just  behind  the  shock. 

The  inviscid  axisymmetric  flow  between  the  cone  surface  and  the  conical  shock  depends 
only  on  the  polar  angle  9.  The  jump  conditions  at  the  shock  must  be  considered  and  the 
velocity  components  may  be  obtained  from  a  numerical  solution  of  the  Taylor-Maccoll  equa¬ 
tion  [8],  Since,  for  a  hypersonic  flow  over  a  slender  cone,  the  density  does  not  vary  much,  we 
use  the  steady,  constant-density  solution  given  by  [9] ,  which  has  the  advantage  of  analytical 
expressions  for  the  velocity  components  and  pressure.  For  a  slender  cone  and  hypersonic 
speeds  these  approximate  solutions  agree  well  with  the  exact  (numerical)  solutions. 

From  the  hypersonic  small-disturbance  approximation  to  the  basic  flow  given  in  [8],  an 
approximate  value  of  the  shock  angle  may  be  obtained  from  the  expression 


sin(6*s  +  9C) 


sin  9C 


i±i+  1  y/2 

2  Mlsm29c) 


(1) 


where  7  is  the  specific  heat  ratio  and  M ^  denotes  the  Mach  number  in  the  free  stream.  For 
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the  weak-shock  solutions,  relevant  to  the  experiments,  this  gives  excellent  agreement  with 
the  exact  solution.  Note  that  the  term  sin  9C  is  0(1).  Thus,  the  shock  angle  may  be 
calculated,  for  a  fixed  cone  angle  and  free-stream  Mach  number.  A  rough  indication  of  the 
results  is  that  9S  ~  9C.  This  is  consistent  with  our  local  analysis  since  we  take  the  thickness 
of  the  upper  deck  (where  the  shock  is  located)  to  be  comparable  with  the  radius  of  the  cone. 
Thus,  the  effect  of  the  shock  should  not  be  neglected.  The  angle  of  the  cone  and  the  flow 
parameters  enable  the  corresponding  value  of  a  to  be  obtained.  Then,  once  the  shock  angle 
has  been  determined,  the  ratio  a/rs  may  be  obtained  simply  from  geometric  arguments.  We 
find  for  a  slender  cone  that 

5  ~  sin^  (2) 

r.5  tan  6S  +  sin  9C ' 

where  we  have  taken  cos  9C  ss  1.  Since  9S  «  9C  a  rough  approximation  gives  a/rs  ~  0.5. 

This  solution  is  not  valid  close  to  the  surface  of  the  cone  so  a  boundary-layer  solu¬ 
tion  has  to  be  introduced  in  this  region.  The  Reynolds  number  of  the  flow  is  defined  by 
Re  =  p-U-L* / /r_.  Taking  the  angle  of  the  cone  to  be  small  the  governing  equations  in  the 
boundary-layer  region  are  given  in  [4]  in  terms  of  dimensionless  coordinates  ( x ,  r,  0)  and 
the  Mach  number,  M ,  just  behind  the  shock.  Then  L*r  is  the  normal  direction  to  the  cone 
surface,  where  f  —  a  on  the  generator  of  the  cone.  The  corresponding  non-dimensional 
velocities  are  (u,  v,  w)  and  the  non-dimensionalised  pressure  and  density  p  and  p,  respec¬ 
tively.  The  boundary  conditions  are  no-slip  at  the  surface  of  the  cone  (coupled  to  the  porous 
layer)  and  appropriate  conditions  at  the  shock  location.  The  non-dimensional  temperature 
and  viscosity  at  the  surface  of  the  cone  are  taken  to  be  Tw  and  pw ,  respectively.  The  only 
restriction  imposed  on  the  temperature  boundary  condition  is  Tw  3>  1,  which  is  violated 
only  for  situations  involving  strong  cooling  on  the  cone  wall  (see  [4].  Usually  the  wall  tem¬ 
perature  is  taken  to  be  Tw  =  TbTr,  where  Tr  is  the  adiabatic  wall  temperature  given  by 
Tr  =  1  +  y/Rr^p/- Ad2  where  Pr  is  the  Prandtl  number.  Thus,  unless  the  constant  Tb  is  very 
small,  Tw  will  be  of  0((y  —  1  )M2)  for  both  adiabatic  walls  (Tb  =  1)  or  isothermal  walls.  The 
analysis  is  unaffected  by  the  particular  choice  of  temperature- viscosity  law.  The  choice  only 

affects  the  bounds  placed  on  various  parameters  of  the  problem.  Sutherland’s  viscosity  law 
—  1  /2  — 

iPw  (1  +  C)TJ  ,  C  being  a  constant)  is  used  henceforth. 


2.1  Porous  boundary 


We  will  present  results  corresponding  to  porous  surfaces  used  in  the  previous  experimental 
investigations.  In  all  cases  the  porous  layer  admittance  Ay  can  then  be  expressed  in  the  form 

Ay  =  -(0o/%)  tanh(A h),  (3) 


where  0 o  is  the  porosity.  The  porous  layer  parameters  non-dimensionalised  with  respect  to 
the  boundary- layer  displacement  thickness  so  h  =  h* / 5*  and  rp  =  r*/8*.  Z0  and  A  are  the 
characteristic  impedance  and  propagation  constant  of  an  isolated  pore,  respectively.  Fedorov 
et  al.  [3]  give  the  following  expressions  for  the  porous  layer  characteristics: 


\/  Pd/Cd 
My/%, 


and  A 


i  toM 

V% 


\J  PdCdi 


(4) 
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where  uj  is  the  disturbance  frequency.  These  are  functions  of  the  complex  dynamic  density 
Pr>  and  complex  dynamic  compressibility  CE.  The  precise  definitions  of  these  quantities 
depends  on  the  structure  of  the  porous  wall  and  are  given  below  for  the  cases  investigated 
here. 

The  wall  boundary  condition,  in  all  cases,  is  then  given  by 

v  =  A, ,(p-p_),  (5) 

where  p_  =  7 -1M-2. 


2.1.1  Regular  microstructure 


Following  [1]  and  [3]  we  consider  the  porous  layer  on  the  cone  surface  to  be  a  sheet  of  thickness 
h*  perforated  with  cylindrical  blind  holes  of  radius  r*  and  equal  spacing  s*  =  This 

model  takes  into  account  gas  rarefaction  effects.  We  have 


Pd  — 


1-F(B„,C) 


,  Cd  =  1  +  (7  —  1)F(Be,  (^y/Pr), 


F(B„,  C)  = 


G{  C) 


1  -  0.5iU2G(C) 


,  F(BEl(VP~r)  = 


G((VP~r) 


}  (6) 


1-0.5 BE({VPr)2G({VPr)’  ) 


where 

Bu  =  (2a“1  —  l)Kn,  BE  =  [7(2q;^1  —  l)/(7  +  l)Pr]Kn,  G(()  =  (7) 


where  (  =  rp\/iupwR/ pw.  Here  J0.i  are  Bessel  functions  of  the  first  kind,  av  and  aE  are 
molecular  accommodation  coefficients,  Kn  is  the  Knudsen  number  and  R  is  the  Reynolds 
number  based  on  boundary-layer  displacement  thickness  of  the  gas  flow. 


2.1.2  Mesh  microstructure 


Following  [6]  we  consider  the  porous  coating  on  the  cone  surface  to  comprise  of  several  layers 
of  stainless  steel  wire  mesh  as  shown  in  figure  3  of  their  paper.  A  similar  model  to  the  one 
described  in  §2.1.1  for  a  regular  microstructure  is  employed.  Following  Kozlov  et  al.  [10]  we 
have  different  expressions  for  the  complex  dynamic  density  and  compressibility.  Following 
[10]  we  can  obtain  the  following  expressions  for  the  porous  layer  characteristics  for  a  square 
mesh  microstructure: 


pD  =  1/(1  -  F( 0),  CD  =  1  +  (7  -  l)F(C),  1 


FIX) 
F(  0 


i+c2E 

m= 0 
00 

i+?£ 

m= 0 


2  /  tanh(/3m)  \ 

llPl  V  An  > 


tanh(/3T, 


> 


(8) 
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where 


(9) 


7m  =  H{m  +  j),  /3m  =  \/(7m“C2),  Pm  =  \J <7m  “  C2)' 

The  characteristic  size  of  an  isolated  pore  is  given  by 

C  =  and  C  =  /PrC-  (10) 

V  dw 

Here  Pr  is  the  Prandtl  number  taken  as  0.72,  a  is  the  pore  width,  R  is  the  Reynolds  number 
based  on  <5*  and  7  is  the  ratio  of  specific  heats  of  air.  The  flow  parameters  are  chosen  to  fit 
the  experimental  conditions  of  [6] . 


2.1.3  Random  microstructure 


Following  [2]  we  consider  the  porous  layer  on  the  cone  surface  to  have  a  random  microstruc¬ 
ture.  A  similar  model  to  the  one  used  for  the  regular  microstructure  is  employed.  We  have 
different  expressions  for  the  complex  dynamic  density  and  compressibility.  Fedorov  et  al.  [2] 
give  the  following  expressions  for  the  porous  layer  characteristics  for  flow  over  a  felt-metal 
microstructure: 


pD  ^00  ( 1  H- 


fl(Al 

Ai 


fi'(A)  —  4/1  + 


daoo/i*  A 


CD  =  7 


7-1 


1  + 


9(^2) 

^2 


cr*(j)0r*2  ’ 

where  the  characteristic  pore  size  is 


=  \2  =  4.Pt\u 

(boa* 


Tid* 


r„  = 


”  2(1  —  </>o)(2  —  <p0) 


(11) 


(12) 


Here  d*  is  the  hbre  diameter,  a*  is  the  flow  resistivity  and  its  value  is  chosen  to  fit  the 
experimental  data  for  flow  over  the  felt  metal.  The  tortuosity  is  taken  to  be  unity. 
Following  Fedorov  et  al.  [2]  rarefaction  effects  are  neglected. 


3  Linear  stability  problem 

The  linear  stability  of  the  basic  flow  described  above  for  M  3>  1  and  Re  3>  1  is  investigated 
in  the  weak-interaction  region  following  a  triple-deck  formulation  used  by  [11]  and  [12].  The 
conditions  to  be  satisfied  at  the  shock  by  a  disturbance  to  this  basic  flow  must  be  specified 
and  these  have  been  derived  in  detail  by  [13].  The  requisite  constraints  were  obtained  by 
considering  the  linearised  jump  conditions  at  the  shock  for  infinitesimal  waves  beneath  the 
shock;  a  similar  procedure  was  adopted  by  [11]  for  flow  over  a  wedge.  Although  the  basic  flow 
is  not  uniform  in  the  regions  below  and  above  the  shock,  [13]  showed  that  the  jump  conditions 
may  still  be  evaluated  at  the  undisturbed  position  of  the  shock.  The  condition  satisfied  by 
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the  pressure  amplitudes  of  the  two  acoustic  waves  (which  are  incident  and  reflected  from  the 
shock)  is  found  to  be  similar  to  that  for  a  wedge  obtained  by  [11]. 

Attention  is  focused  at  a  location  on  the  surface  of  the  cone  with  non-dimensional  radius 

3  i  _ :£  _ 9 

a  =  a* /L*.  It  is  assumed  that  aRe* M 3 fj,w 8 Tw  8  =  a  ~  0(1)  denotes  the  scale  of  the  radius 

_ 3  _i  3  9 

at  this  point;  thus  we  have  chosen  sindc  ~  9C  ~  Re  8M  ifiwTw.  We  note  that  in  terms 
of  the  interaction  parameter  y  =  M^a* / L*  ~  0(1)  corresponding  to  a  moderate  inviscid 
interaction.  Since  from  (1),  9S  ~  9C,  the  shock  is  located  in  the  upper  deck  of  the  triple-deck 
structure. 


Our  study  is  confined  to  the  question  of  the  stability  of  the  flow  at  a  location  on  the 
body  where  the  boundary- layer  thickness  is  0(i?e_5 L*),  which  is  thin  compared  to  the  local 
radius  of  the  cone.  This  situation  is  chosen  so  that  curvature  effects  are  significant.  The 
analysis  is  somewhat  simplified  if  non-parallel  effects  can  be  neglected  and  [11]  showed  that 
this  is  justifiable  if  the  ‘Newtonian’  assumption  7  —  1  -C  1  is  made.  Thus,  for  simplicity, 
this  condition  is  taken  to  hold  in  the  following  analysis  although  it  can  be  easily  relaxed  for 
more  involved  studies. 

It  is  convenient  to  scale  out  some  of  the  parameters  in  the  problem,  namely  fiw ,  Tw  and 
A,  where  the  last  quantity  denotes  the  boundary-layer  skin  friction.  For  axisymmetric  flow 
the  Mach  number  may  be  scaled  out  of  the  linear  stability  problem. 

We  consider  perturbations  proportional  to 


E  =  exp  [i  ( aX  +  ncf)  —  fir)] , 


where  a  and  n  are  the  streamwise  and  azimuthal  wavenumbers  respectively  and  fl  is  the 
frequency  of  the  disturbance.  Note  that  n  is  an  integer. 

Previous  scaling,  following  [11]  and  [14]  is  applied  to  the  resulting  equations  for  axisym¬ 
metric  disturbances  (corresponding  to  n  —  0).  Analytic  solutions  of  these  equations  yields 
an  eigenrclation  relating  the  streamwise  wavenumber  a  and  frequency  fl,  namely 

Ai'QCo)  _  \  1/3  rj-  .  \  Io(iars)K0(iaa)  -  I0(iaa)K0(iars) 

/^Ai(£)d£  5  I0(iars)Kl(iaa)  +  Ii(iaa)K0(iars) 

Here  £0  =  — i ^flar2/3,  Ai(£)  is  the  Airy  function,  Kn(z )  and  In{z)  are  the  usual  modified 
Bessel  functions,  and  Ay  =  Re~1^s/j,w8\1^4Tw8(M2  —  1  )3//8kLy.  The  admittance,  Ay,  is  a 
function  of  the  disturbance  frequency  and  depends  on  the  physical  properties  of  the  flow 
and  the  porous  layer.  The  angular  frequency  of  disturbance  propagation  through  the  pore 
is  u  =  (R/Re)Re1/4Hw1^4\3/2Tw3/4(M2  —  1)  1//40,  where  R  is  the  Reynolds  number  based 
on  displacement  thickness.  The  parameter  rs  is  the  scaled  non-dimensional  location  of  the 
shock.  The  values  of  a  and  rs  depend  on  the  physical  parameters  for  the  flow. 


A  similar  analysis  is  applied  for  non-axisymmetric  disturbances  with  azimuthal  wavenum¬ 
ber  n.  The  scaling  is  different  and  in  particular  Ay  =  Re-1/8 /mJ8 A1//4T^8 M5/4Ay  and 

ui  =  (R/Re)Re1/4nw1/4\3/2Tw3/4M-1/2Q.  The  resulting  eigenrclation  for  non-axisymmetric 
disturbances  is  given  by 


Ai'(6) 

JTAim 


(k*)1/3 


Ay  + 


In(iars)Kn(iaa ) 
In(iars)K!n{iaa) 


In(iaa)Kn(mrs) 
I'n(\aa)Kn(\ars) ' 


(14) 
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In  order  to  present  our  results  in  regimes  of  practical  interest,  we  will  use  the  flow 
parameters  from  the  relevant  experimental  studies  [1,2]  and  [6].  The  cone  angle  and  Mach 
number  from  the  experiments  will  determine  the  shock  angle  6S  and  the  scaled  radius  a.  Then 
the  scaled  shock  location  rs  is  determined  as  described  above.  The  experimental  conditions 
correspond  to  a/rs  =  0.57  thus,  we  choose  to  present  our  results  for  this  case. 

The  values  for  flow  of  a  perfect  gas  were  chosen  as  M  =  5.3,  T.w  =  Tad,  7  =  1.4,  Pr  =  0.71, 
Re  1  =  15.2  x  106,  T*  =  56.4 K  and  au  =  cxe  =  0.9.  These  values  correspond  to  the  ones 
used  in  the  previous  numerical  studies  of  [15].  The  experiments  were  conducted  on  a  cone 
of  0.5m  in  length. 


4  Neutral  results 

We  now  consider  neutrally  stable  solutions  of  the  dispersion  relations  (13)  and  (14)  corre¬ 
sponding  to  the  axisymmetric  case  and  the  non-axisymmetric  case,  respectively.  We  seek 
solutions  with  a  and  fl  real  and  investigate  the  effect  of  the  porous  microstructure  on  the 
neutral  solutions.  The  presence  of  the  shock  allows  for  multiple  modes  of  solutions. 

4.1  Regular  microstructure 

The  current  results  for  the  regular  microstructure  comprising  a  regular  array  of  cylindrical 
pores  of  circular  cross-section  are  compared  to  the  results  for  a  solid  wall  for  axisymmetric 
and  non-axisymmetric  modes. 

For  the  regular  microstructure  of  [3]  the  results  presented  below  are  for  the  porous  layer 
parameters  r*  =  28.5 gm,  (f>0  =  0.2  and  h*  »  r*.  The  last  relation  implies  that  Ah  — >  00 
so  Eq.  (3)  simplifies  to  Ay  =  —(j)0/Z0.  The  effect  of  pore  depth  is  investigated  in  section  6. 
The  porous  layer  characteristics  correspond  to  experimental  values  used  in  [3]. 

The  flow  conditions  match  the  experimental  conditions  of  [6] .  Here  a  cone  of  angle  9C  =  7° 
was  used  with  a  flow  of  Mach  number  =  6. 

The  effects  of  rarefaction  are  included  and  the  Knudsen  number  is  calculated  from  the 
expression 

Kn  =  -  yfin ryTw, 

r  PR 

giving  Kn  =  0.494. 

Figure  2  shows  the  neutral  values  of  frequency  H  as  a  function  of  radius  a  for  the  Erst 
four  axisymmetric  modes  for  a/rs  =  0.57.  The  scaled  radius  varies  linearly  with  L*.  The 
range  0  <  a  <  5  corresponds  to  0  <  L*  <  0.02m  for  axisymmetric  modes  for  all  the  wall 
models  considered.  The  dashed  lines  are  for  the  porous  wall  and  the  solid  lines  correspond 
to  the  results  for  a  solid  wall  from  [4],  The  flow  is  unstable  in  regions  above  these  curves. 
For  the  higher  modes  for  the  solutions  for  a  there  is  no  discernible  difference  between  the 
results  (not  shown).  However,  the  neutral  solutions  for  are  lower  than  those  for  a  solid 
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Figure  2:  Neutral  values  of  frequency  for  the  first  four  axisymmetric  modes  as  a  function 
of  radius  a  with  n= 0  and  a/rs  =  0.57: - regular  porous  wall; - solid  wall. 

wall,  particularly  for  larger  values  of  a.  Thus,  the  flow  over  the  porous  surface  will  become 
unstable  for  lower  frequencies  than  those  for  the  solid  wall. 

Figure  3  shows  the  neutral  solutions  for  hi  as  a  function  of  a  for  non-axisymmetric  dis¬ 
turbances  for  a/rs  =  0.57  and  n  —  1  and  n  —  2  for  a  solid  wall  and  a  regular  microstructure. 
For  n  —  1  the  range  0  <  a  <  5  corresponds  to  0  <  L*  <  14m.  Again  the  porous  wall  has 
only  a  small  effect  on  the  neutral  values  of  a  (not  shown).  The  behaviour  of  the  first  mode 
is  different  from  the  higher  modes  with  f!  — >  0  as  a  — »  0.  We  see  from  figure  3  that  the 
porous  wall  has  a  very  small  effect  on  the  neutral  values  of  for  this  first  mode.  However, 
the  porous  wall  has  a  significant  effect  on  the  higher  modes,  leading  to  lower  neutral  values 
of  fh  Thus,  the  effect  of  the  porous  wall  is  destabilising.  Moreover,  this  effect  increases  for 
the  higher  modes.  Similar  results  have  been  obtained  for  n  =  3. 

4.2  Mesh  microstructure  and  random  microstructure 

The  neutral  results  for  the  felt-metal  microstructure  are  compared  to  those  for  a  mesh 
microstructure.  Following  [6]  we  consider  the  mesh  microstructure  to  correspond  to  0o  =  0.8 
and  a  =  0.05mm.  In  the  experiments  the  porosity  of  the  mesh  microstructure  is  larger  than 
that  of  the  regular  array  of  circular  pores.  This  will  lead  to  a  larger  destabilising  effect  on 
the  first  Mack  mode  (see  section  5).  We  show  our  results  for  h*  a. 

Figure  4  shows  the  neutral  values  of  H  as  a  function  of  a  for  the  mesh  microstructure 
and  the  felt-metal  microstructure  for  the  first  four  axisymmetric  modes.  We  find  that  the 
felt  metal  has  a  larger  destabilising  effect  than  the  mesh  microstructure,  resulting  in  lower 
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Figure  3:  Neutral  values  of  frequency  for  the  first  five  non-axisymmetric  modes  as  a 

function  of  radius  a  with  a/rs  =  0.57: - regular  porous  wall; - solid  wall;  (a)  n  =  1, 

(b)  n  =  2. 


values  of  neutral  frequency. 

Figure  5  shows  the  neutral  solutions  for  as  a  function  of  a  for  non-axisymmetric  dis¬ 
turbances  for  n  —  1  and  a/rs  =  0.57  for  a  mesh  microstructure  and  a  random  microstructure 
for  the  first  five  modes.  The  neutral  values  of  f2  for  the  mesh  microstructure  are  lower  than 
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Figure  4:  Neutral  values  of  frequency  for  the  first  four  axisymmetric  modes  as  a  function  of 
radius  a  with  n— 0  and  a/rs  =  0.57: - mesh  microstructure; - random  microstructure. 


Figure  5:  Neutral  values  of  frequency  for  the  first  five  non-axisymmetric  modes  as  a 

function  of  radius  a  with  n—  1  and  a/rs  =  0.57: - mesh  microstructure;  - random 

microstructure. 
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those  for  the  regular  microstructure  comprising  circular  pores  as  a  result  of  the  porosity  being 
larger.  We  see  from  figure  5  that  the  random  microstructure  has  a  significant  destabilising 
effect,  with  the  neutral  values  much  lower  than  those  for  the  mesh  microstructure. 


5  Spatial  Growth  Rates 

We  now  concern  ourselves  with  an  examination  of  the  spatial  evolution  of  disturbances  so 
that  we  concentrate  on  solutions  of  Eq.  (13)  and  Eq.  (14)  with  fl  real  and  a  complex.  If 
ol  =  ar  +  icq  then  cq  >  0  is  indicative  of  stability  while  cq  <  0  denotes  spatial  instability. 

Figures  6-8  show  the  dependence  of  the  spatial  growth-rate  parameter  cq  on  the  mode 
frequency  for  a  regular  microstructure  compared  to  a  solid  wall  for  n  —  0,  n  —  1  and  n  =  2, 
respectively.  The  results  are  shown  for  a/rs  =  0.57  and  a  =  0.6,  1.0,  1.5  and  2.0.  There  is  a 
complete  family  of  modes,  as  we  saw  in  onr  account  of  neutral  disturbances,  and  it  is  clear 
that  for  each  member  of  the  family  there  is  a  cut-off  frequency  flc  such  that  for  that 

particular  mode  is  stable  but  it  becomes  unstable  if  O  >  fic.  The  effect  of  the  neutral  values 
being  altered  by  the  presence  of  the  porous  wall  for  all  wall  models  leads  to  variation  in  the 
growth  rates.  The  maximum  values  of  (— cq)  for  the  porous  wall  are  increased  greatly  from 
those  of  a  solid  wall  for  each  mode. 

In  addition,  we  also  observe  that  the  growth  rate  for  the  porous  modes  does  not  rapidly 
approach  zero  at  high  frequencies,  as  is  the  case  for  the  solid  wall.  Figure  6  shows  that  for 
the  porous  wall  and  the  solid  wall  the  first  mode  has  the  largest  growth  rate  for  axisymmetric 
disturbances  for  smaller  values  of  a.  However,  for  larger  values  of  a  the  second  mode  has  the 
largest  maximum  growth  rate.  As  a  increases  the  growth  rates  decrease. 

From  figure  7  we  see  that  for  non-axisymmetric  modes  with  n  —  1  the  growth  rate  of 
the  higher  modes  is  greatly  enhanced  by  the  regular  porous  wall.  It  is  the  last  mode  shown 
which  has  the  highest  growth  rate.  The  results  for  cq  presented  in  figure  8  for  n  —  2  are 
similar  to  those  for  n  —  1. 

Figure  9  shows  the  results  for  cq  for  the  felt-metal  microstructure  with  0o  =  0.75.  These 
are  compared  to  the  values  for  a  mesh  microstructure  with  0o  =  0.8.  The  results  are 
presented  for  n  =  0,  a/rs  =  0.57  and  a  =  0.6,  1.0,  1.5,  2.0.  We  find  that  the  growth  rates 
of  all  modes  are  increased  greatly  for  the  felt-metal  microstructure  compared  to  the  mesh 
microstructure.  For  the  felt-metal  microstructure  the  first  mode  has  the  largest  growth  rate, 
switching  to  the  second  mode  as  the  radius  a  increases.  This  occurs  at  a  larger  value  of  a 
than  for  the  mesh  microstructure.  These  results  agree  with  numerical  results  for  a  random 
microstructure  from  [7]. 

Figure  10  shows  the  spatial  growth  rates  for  n  —  1  for  the  mesh  microstructure  and  the 
felt-metal  microstructure.  The  effect  of  the  felt-metal  microstructure  is  to  give  a  very  large 
increase  in  the  maximum  growth  rates,  approximately  four  times  in  magnitude,  compared 
to  the  mesh  microstructure. 

We  end  this  discussion  on  the  effects  of  a  porous  wall  on  the  spatial  growth  rates  by 
considering  the  effects  of  rarefaction  for  a  regular  microstructure  comprising  an  array  of 
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Figure  6:  Spatial  growth-rate  parameters  ctj(fl)  for  the  first  few  non- neutral  axisymmetric 

modes  for  a/rs  =  0.57:  (a)  a  =  0.6;  (b)  a  =  1.0;  (c)  a  =  1.5;  (d)  a  =  2.0  : - regular 

porous  wall; - solid  wall. 


circular  pores.  In  figure  11  we  show  the  neutral  values  of  11  as  a  function  of  a  for  the  first 
five  non-axisymmetric  modes  with  n  =  1,  a/rs  =  0.57  for  Kn  =  0  and  Kn  =  0.494.  We  see 
that  the  effect  of  non-zero  Knudson  number  is  destabilising,  with  lower  values  of  11  compared 
to  Kn  =  0,  particularly  for  the  higher  modes.  The  effect  of  rarefaction  on  the  spatial  growth 
rates  is  illustrated  in  figure  12,  showing  the  results  corresponding  to  figure  11.  The  growth 
rates  are  increased  for  non-zero  Knudson  number,  with  larger  increases  for  the  higher  modes. 
We  note  that  for  the  Mack  second  mode  [8]  found  that  the  effect  of  rarefaction  is  stabilising. 


6  Maximum  Spatial  Growth  Rates 

We  investigate  the  effects  of  changing  the  properties  of  the  porous  wall  on  the  linear  stability 
of  hypersonic  flow  over  a  sharp  cone.  We  choose  to  concern  ourselves  with  an  examination 
of  the  effect  of  the  physical  properties  of  the  porous  layer  on  the  maximum  spatial  evolution 
of  disturbances.  Thus,  we  consider  the  mode  which  has  the  largest  value  of  (—«;). 

The  effects  of  porosity,  pore  radius  and  pore  depth  on  the  maximum  spatial  growth  rates 
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Figure  7:  Spatial  growth-rate  parameters  aj(fl)  for  the  first  few  non- neutral  non- 
axisymmetric  modes  for  a/rs  =  0.57  and  n  =  1:  (a)  a  =  0.6;  (b)  a  =  1.0;  (c)  a  =  1.5; 
(d)  a  =  2.0  : - regular  porous  wall; - solid  wall. 

are  considered  for  azimuthal  wavenumber  n  —  1  and  a/rs  =  0.57  and  a  =  0.8. 

The  effect  of  increasing  the  porosity  0o  (by  decreasing  the  pore  spacing)  for  fixed  pore 
radius  r*  =  30 firm  is  illustrated  in  figure  13.  We  see  that  as  the  porosity  increases  the 
maximum  spatial  growth  rate  amax  increases.  Here  Umax  is  the  maximum  value  of  —  for 
the  different  modes  for  particular  values  of  a/rs. 

In  figure  14  we  show  the  effect  of  the  pore  radius,  r* ,  on  the  maximum  spatial  growth 
rate  Umax  for  0o  =  0.25.  We  see  that  larger  pore  radius  leads  to  larger  maximum  spatial 
growth  rates. 

The  previous  results  presented  have  been  obtained  by  assuming  an  infinite  pore  depth. 
The  effect  of  finite  pore  depth  is  demonstrated  in  figure  15  for  porosity  0o  =  0.25  and 
r*  =  30 pm.  We  see  that  as  h/rp  increases  the  maximum  spatial  growth  rate  increases  but 
it  levels  off  for  h/rp  >  20.  In  the  experiments  the  perforated  sheet  was  relatively  thick 
with  h*  =  0.45-0. 5mm  and  r*  =  25-30/jm,  giving  h/rp  =  15-20.  This  verifies  our  earlier 
assumption  of  taking  h/rp  3>  1  in  our  calculations.  Previous  calculations  determined  that 
maximum  stabilisation  of  the  second  Mack  mode  was  achieved  for  h/rp  =  3. 
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Figure  8:  Spatial  growth-rate  parameters  a*  (12)  for  the  first  few  non- neutral  non- 
axisymmetric  modes  for  a/rs  =  0.57  and  n  =  2:  (a)  a  =  0.6;  (b)  a  =  1.0;  (c)  a  =  1.5; 
(d)  a  =  2.0  : - regular  porous  wall; - solid  wall. 

In  figure  16  we  show  the  effect  of  pore  depth  on  a  mesh  microstructure.  An  increase 
in  the  depth  corresponds  to  more  layers  of  the  fine  mesh.  We  find  that  for  h*  =  0.15mm 
(corresponding  to  three  layers  of  mesh)  the  maximum  growth  rate  is  a  maximum  for  the 
first  Mack  mode.  The  previous  results  of  [6]  showed  that  h*  =  0.15mm  gave  the  maximum 
stabilisation  of  the  second  Mack  mode.  So  this  is  the  worst  case  for  the  first  Mack  mode. 

In  figure  17  we  present  the  maximum  spatial  growth  rates  to  correspond  to  the  exper¬ 
imental  conditions  of  Maslov  [15]  for  a  regular  microstructure  with  r*  =  25/mi,  (j)0  =  0.2 
and  Re  =  10  x  106.  Here  the  dimensional  spatial  growth  rate  is  shown  as  a  function  of  the 
distance  along  the  cone  for  the  mode  giving  the  largest  growth  rate  for  n  =  0,  1,  2,  3  and 
a/rs  =  0.57.  The  dashed  lines  correspond  to  the  porous  wall  and  the  solid  lines  correspond  to 
a  solid  wall.  We  see  that  for  the  non-axisymmetric  modes  the  growth  rates  are  significantly 
larger  for  the  porous  wall  compared  to  the  solid  wall.  However,  the  size  of  the  growth  rates 
is  much  smaller  than  those  obtained  for  the  second  Mack  mode  in  [15]. 
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Figure  9:  Spatial  growth-rate  parameters  for  the  first  few  non- neutral  axisymmetric 

modes  for  a/rs  =  0.57:  (a)  a  =  0.6;  (b)  a  =  1.0;  (c)  a  =  1.5;  (d)  a  =  2.0: - mesh 

microstructure; - random  microstructure. 

7  Nonlinear  analysis 

The  linear  stability  analysis  will  not  be  valid  for  larger  disturbances.  Thus,  it  is  important 
to  determine  the  effect  of  nonlinearity  on  the  stability  of  hypersonic  boundary  layer  flow 
over  a  sharp  slender  cone  with  a  porous  wall.  This  has  been  investigated  experimentally 
for  second  mode  disturbances  using  the  bicoherence  method  by  Chokani  et  al.  [16]  and 
Bountin  et  al.  [17],  We  extend  the  linear  analysis  into  the  weakly  nonlinear  regime.  Here 
larger  disturbances  are  considered  but  not  so  large  as  to  alter  the  basic  flow.  For  the  current 
problem  with  a  solid  wall  this  has  been  investigated  by  Stephen  [5],  The  planar  case  was 
considered  previously  by  Seddougui  and  Bassom  [18].  The  lengthy  analysis  for  the  hypersonic 
flow  over  a  cone  with  a  porous  closely  follows  these  studies.  Thus,  the  main  results  will  be 
summarised  and  the  differences  highlighted. 

We  consider  a  weakly  nonlinear  disturbance  which  develops  in  the  vicinity  of  a  linear 
neutral  point  (real  a  and  f!  for  fixed  n  >  0).  If  the  relative  amplitude  of  the  disturbance  in 
the  lower  deck  of  the  triple-deck  structure  is  O(h),  h  <C  1,  then  the  scaled  amplitude  A  of 
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Figure  10:  Spatial  growth-rate  parameters  cq(fi)  for  the  first  few  non- neutral  non- 
axisymmetric  modes  for  a/rs  =  0.57  and  n  =  1:  (a)  a  =  0.6;  (b)  a  =  1.0;  (c)  a  =  1.5; 
(d)  a  =  2.0: - mesh  microstructure; - random  microstructure. 

the  mode  will  evolve  on  an  0(h2)  lengthscale.  Thus,  we  consider  perturbations  at  the  point 

x  —  1  +  h2x  2- 

Since  the  skin  friction  A  is  a  function  of  x,  it  will  also  be  slightly  perturbed  from  its  neutral 
value  and  we  write 

A  =  1  +  h~\2, 

where  A2  =  x2d\/ dx\x=\.  Additionally,  we  choose  to  fix  the  azimuthal  wavenumber  and 
write 

O  =  -F  Ji2 fl2, 

where  Hi  is  the  neutral  value  of  the  frequency  obtained  from  the  linear  stability  problem.  To 
account  for  the  slow  modulation  of  the  amplitude  on  streamwise  lengthscales  we  introduce 
the  coordinate 

A  =  h2X, 

and  then  by  use  of  multiple  scales  we  replace  all  X  derivatives  throughout  according  to 

d  d  2  d 
dX  dX  +  1  W 
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Figure  11:  Neutral  values  of  frequency  fl  for  the  first  five  non-axisymmetric  modes  as  a 
function  of  radius  a  for  a  regular  microstructure  with  n—  1  and  a/rs  =  0.57: - Kn=0.494; 


Kn=0. 


Figure  12:  Spatial  growth-rate  parameters  cq(fl)  for  the  first  few  non- neutral  non- 
axisymmetric  modes  for  a  regular  microstructure  for  a/rs  =  0.57,  n  —  1  and  a  =  0.6:  - 
—  Kn=0.494; - Kn=0. 
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Figure  13:  Maximum  spatial  growth  rates  Umax  for  varying  porosity  (j) o  for  a  regular  mi¬ 
crostructure  for  n  =  1,  a/rs  =  0.57,  r*  =  30 /im  and  a  =  0.8. 


Figure  14:  Maximum  spatial  growth  rates  Umax  f°r  varying  pore  radius  r*  for  a  regular 
microstructure  for  n  =  1,  a/rs  =  0.57,  0 o  =  0.25  and  a  =  0.8. 


Now  for  /*  <C  1  we  proceed  to  seek  solutions  of  the  lower-deck  equations  ((8)  and  (9)  in  the 
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Figure  15:  Maximum  spatial  growth  rates  Umax  for  varying  pore  depth  ratio  h/rp  for  a 
regular  microstructure  for  n  =  1,  a/rs  =  0.57,  0o  =  0.25,  r*  =  30pm  and  a  =  0.8. 


Figure  16:  Maximum  spatial  growth  rates  <rmax  for  varying  pore  depth  h*  for  n  =  1, 
a/rs  =  0.57,  (f>0  =  0.25  and  a  =  0.8  for  a  mesh  microstructure  for  a  =  0.05 mm. 

second  report)  and  the  upper-deck  problem  ((11)  in  the  second  report)  where 

3 

(U-(l  +  h2 A2)  y,  fo  W,  P,  A,p)=J2  hj  (Uj,  Vj7  Wj,  Aj/pj)  +  0(h4).  (15) 

3= 1 
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Figure  17:  Maximum  spatial  growth  rates  Umax  as  a  function  of  L*  for  a  regular  microstruc¬ 
ture  for  n  =  0,  1,  2,  3,  a/rs  =  0.57,  0o  =  0.2  and  r*  =  25 /im. 

Substitution  of  (15)  into  the  disturbance  equations  leads  to  a  hierarchy  of  problems  at 
increasing  orders  in  h. 

The  solution  for  the  O(h)  terms  in  (15)  has  been  obtained  from  the  linear  stability 
analysis  resulting  in  the  dispersion  relations  (13)  and  (14).  The  solution  for  A\  has  the  form 

A\  =  AnE  +  A\Ie  1, 

where  the  superscript  (c)  denotes  complex  conjugate. 

7.1  The  second-order  problem 

At  0(h2)  we  find  the  solution  for  A2  takes  the  form 

A2  =  A22E-  +  A2q  +  A22  E  2, 

with  similar  expansions  for  U2)  V2,  W2,  P2  and  p2.  The  analysis  follows  that  for  the  solid  wall 
problem  (see  Smith  [19])  but  the  boundary  condition  at  the  porous  wall  leads  to  additional 
terms  in  the  solutions.  Solving  the  equations  and  satisfying  all  the  boundary  conditions  leads 
to  an  expression  relating  the  amplitude  at  0(h 2)  to  the  square  of  the  amplitude  at  0(h). 

7.2  The  third-order  problem 

At  this  order  the  amplitude  equation  for  the  unknown  function  An(A^)  is  determined.  This 
equation  arises  from  a  solvability  condition  on  the  0(h3E)  terms.  We  seek  solutions  of  the 
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form 


p3  =  p31E  +  p32E 2  +  p33E3  +  p$E  1  +  p$E  2  +  p33  E  3 +  p30. 

The  solution  in  the  lower  deck  involves  an  inhomogeneous  differential  equation  of  which  a 
solution  only  exists  if  a  certain  compatibility  condition  holds.  This  condition  is  derived  by 
considering  the  adjoint  system  of  the  problem  as  in  Hall  and  Smith  [20],  Specifically,  we 
multiply  the  inhomogeneous  equation  by  the  adjoint  function  and  integrate  over  the  range 
[fo,oo]. 

By  matching  the  solution  for  p33  with  the  solution  for  the  disturbance  pressure  in  the 
main  deck  we  obtain  the  evolution  equation  for  An  in  the  form 


«i  ~  —  (02-^2  +  03^2)  An  +  a4H11|H1i|2 
dA 


The  complex  constants  in  (16)  are  given  by 


1  "2  Ai'(&) 

-  i 


K  Otcr  K 


Pi  /  [^(OT(O]^-2M1/3Ai(e0)aPr1 

•J  £0 


(K  fn  +  Cndn)  ( Kn(iaa)  -  Uaj  1/3P iK^iaa) 

—{bndn  +  cnen)  ( In(iaa )  -  A‘  —  (io')~1/3Pi J^(io-a) 


a2  — 

i  2^3aL2  +  cm  1 

a3  = 

(m*)1/3k_1T18, 

r2/3  (a)5/3 

r 

a>4  = 

1  1 0 

All 

k\k\z 

- 

(16) 


(17) 

(18) 
(19) 


ill  -  Ay  [a2/3PiT22  -  P2T25L12]  -  ylv'^pfiisl  ,  (20) 


where  k  =  f™Ai(s)ds,  Pi  =  (Hy  +  ^)-1,  P2  =  (Hy  +  2^)-1.  All  other  abbreviations  Tij, 
Lij  are  defined  in  Appendix  A.  The  coefficients  bn ,  cn,  dn,  en,  fn  depend  on  the  neutral  values 
of  a  for  fixed  values  of  a  and  rs  and  involve  modified  Bessel  functions.  They  are  defined 
in  the  appendix  of  Stephen  [5].  The  remaining  complex  constants  involve  Airy  functions 
and  are  defined  in  Appendix  A.  The  coefficients  are  functions  of  the  neutral  results  and  so 
depend  on  the  porous  wall  properties.  The  porous  wall  has  introduced  additional  terms  in 
the  coefficients.  Thus,  their  evaluation  requires  substantial  numerical  calculations. 

The  corresponding  nonlinear  problem  for  axisymmetric  disturbances  must  be  considered 
separately  since  the  Mach  number  can  be  scaled  out  of  the  weakly  nonlinear  problem.  The 
analysis  is  very  similar  to  that  for  the  non-axisymmetric  problem  and  the  resulting  amplitude 
equation  is  presented  in  Appendix  B. 


8  Nonlinear  results 

We  will  begin  this  section  by  presenting  the  results  of  the  nonlinear  stability  analysis  for 
the  solid  wall  with  flow  parameters  corresponding  to  the  experiments  of  [15]  and  [3].  This 
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was  not  done  in  [5].  We  will  then  present  corresponding  results  obtained  using  the  different 
porous  wall  models  and  deduce  the  effect  of  the  porous  walls  on  the  nonlinear  stability. 


8.1  Solid  wall 


In  order  to  present  our  results  in  regimes  of  practical  interest,  we  will  use  the  flow  parameters 
from  the  relevant  experimental  studies.  Thus  all  our  results  are  presented  for  a/rs  =  0.57. 

We  turn  to  the  non-axisymmetric  problem  and  the  evolution  equation  (16).  We  can  solve 
this  equation  using  separation  of  variables.  Following  [21]  we  can  determine  an  explicit 
expression  for  the  amplitude  |An|2  as 


|Ai|2 


2  K,e™ 
KCx  -  2Re 


where  /C  =  2Re  (  — A2  +  — ^2  )  and  C\  is  a  constant  of  integration.  When  JC  >  0  we  have 
linear  instability.  Note  that  A2  is  negative  downstream  of  the  neutral  location.  If  then 
Re  <  0,  nonlinear  effects  are  stabilising  and  the  linearly  unstable  mode  is  supercritically 
stable  with  an  equilibrium  amplitude  given  by 


|Alii 


Re  (S^  + 


\ 


In  figure  18a  we  show  Re(a2/oq)  as  a  function  of  a  for  n  —  1  corresponding  to  the  first  five 
neutral  modes  of  (14)  for  a  solid  wall.  The  arrows  here  and  on  all  the  subsequent  figures 
indicate  increasing  mode  number.  We  can  see  that  this  quantity  is  always  negative.  The 
results  for  n  —  2  are  shown  in  figure  19a.  The  magnitudes  are  decreased  when  compared  to 
n  —  1 .  The  corresponding  values  of  Re(as/ai)  are  shown  in  figures  18b  and  19b,  respectively 
for  n  =  1  and  n  =  2.  There  is  a  difference  in  behaviour  of  Re(a3/ai)  for  the  first  mode  for 
small  values  of  a,  corresponding  to  the  anomalous  behaviour  of  the  lowest  neutral  solution; 
see  figure  3.  The  effect  of  increasing  the  azimuthal  wavenumber  is  to  increase  the  magnitude 
of  Re(a3/ai). 

We  now  investigate  the  effect  of  nonlinearity  on  linearly  unstable  disturbances  by  consid¬ 
ering  the  sign  of  R,e(a4/oi).  In  order  to  see  the  behaviour  of  the  different  modes  the  value  of 
Re(a4/ai)  versus  a  for  0  <  a  <  1  is  shown  in  figure  20a  and  versus  a  for  1  <  a  <  5  is  shown  in 
figure  20b  for  n  =  1.  In  figures  20a  and  20b  we  can  see  that  the  sign  of  Re(a4/ai)  is  always 
positive  for  the  first  mode  (m  =  1).  Thus  nonlinear  effects  always  destabilise  this  mode 
possibly  leading  to  a  subcritical  instability.  The  effect  of  nonlinearity  on  the  remaining  four 
modes  depends  on  the  value  of  a.  For  a  <  1,  Re(a4/ai)  <  0  leading  to  supercritical  instabil¬ 
ity.  As  the  value  of  a  increases,  the  sign  of  Re(a4/ai)  becomes  positive  beginning  with  the 
higher  modes  indicating  that  nonlinear  effects  now  destabilise  these  linearly  unstable  modes. 
In  the  limit  of  large  a  we  can  see  that  Re^/eq)  — »•  0  for  all  the  modes,  with  the  first  mode 
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Figure  18:  (a)  Re(a2/oi);  (b)  Re(a3/ai),  as  a  function  of  local  cone  radius  a  for  the  first  five 
modes.  Results  are  shown  for  solid  wall  with  n  =  1  and  a/rs  =  0.57. 


Figure  19:  (a)  Re(a2/ai);  (b)  Re(a3/ai),  as  a  function  of  local  cone  radius  a  for  the  first  five 
modes.  Results  are  shown  for  a  solid  wall  with  n  —  2  and  a/rs  =  0.57. 
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a 


Figure  21:  (a)  Re(a4/ai)  for  the  first  five  modes  as  a  function  of  local  cone  radius  a  for 
0  <  a  <  2.3.  Results  are  shown  for  a  solid  wall  with  n  —  2  and  a/rs  =  0.57;  (b)  Re(a4/ai) 
for  the  hrst  mode  as  a  function  of  a. 


Figure  22:  Re(a4/ai)  for  the  hrst  hve  modes  as  a  function  of  local  cone  radius  a  for  2.3  < 
a  <  5.0.  Results  are  shown  for  a  solid  wall  with  n  —  2  and  a/rs  =  0.57. 
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having  the  highest  amplitude.  For  a  fixed  cone  angle  the  effect  of  increasing  a  is  to  move 
further  along  the  cone  surface.  Thus  at  large  streamwise  distances  we  can  expect  nonlinear 
disturbance  amplitudes  to  be  very  small.  The  corresponding  results  for  Re(a4/a4)  with  n  —  2 
can  be  seen  in  figures  21  and  22.  In  figure  21a  we  can  see  that  effect  of  nonlinearity  on  mode 
numbers  >  1  for  n  —  2  are  similar  to  that  for  n  =  1.  Higher  azimuthal  wavenumber  allows 
the  stabilising  effect  of  nonlinearity  to  persist  for  larger  ranges  of  a.  In  figure  22  we  see  that 
Re(a4/ai)  becomes  positive  for  a  >  2.3  with  the  higher  modes  becoming  destabilised  first. 
The  effect  of  azimuthal  wavenumber  is  more  significant  on  the  first  mode  as  can  be  seen  in 
figure  21b.  Here  we  see  that  the  first  mode  is  stabilised  by  nonlinearity  for  a  narrow  range 
of  a  (Re(a4/ai)  <  0  when  1.6  <  a  <  3.3).  In  the  limit  of  large  a  nonlinear  effects  on  the 
first  mode  persist  while  the  effects  on  the  higher  modes  diminish  as  Re(a4/ai)  — >  0  for  these 
modes.  The  effect  of  nonlinearity  for  disturbances  with  n  —  3  have  also  been  investigated. 
The  results  (not  shown)  indicate  that  overall  trends  remain  similar  to  that  for  n  =  2.  We 
may  thus  conclude  that  nonlinear  effects  tend  to  stabilise  the  higher  modes  for  a  wider  range 
of  a  for  higher  azimuthal  wavenumbers,  and  the  first  mode  becomes  the  most  destabilised 
by  nonlinearity  at  large  values  of  a. 

We  consider  the  effect  of  nonlinearity  on  axisymmetric  disturbances  next.  The  neutral 
linear  stability  results  for  these  disturbances  were  shown  in  figure  2.  Axisymmetric  distur¬ 
bances  are  linearly  unstable  if  Re^o-Waio)  >  0.  If  Re(a40/ai0)  <  0  for  these  disturbances 
then  nonlinear  effects  are  stabilising  and  the  linearly  unstable  modes  are  supercritically  stable 
with  an  equilibrium  amplitude 


\An 


/— Re(q20)\  1/2 
Re(a40)  / 


In  figure  23a  we  can  see  Re(a2o/oio)  as  a  function  of  a.  We  notice  that  this  quantity  is  always 
negative.  Figure  23b  shows  Re(a40/ai0)  versus  a  corresponding  to  the  first  four  modes.  We 
see  that  this  quantity  is  negative  for  all  the  modes  with  the  exception  of  the  first  mode.  For 
this  mode  R,e(«4o/«io)  becomes  slightly  positive  for  a  >  2.3.  Thus,  we  can  expect  nonlinear 
effects  to  stabilise  linearly  unstable  axisymmetric  disturbances  with  the  exception  of  the  first 
mode  which  is  slightly  destabilised  above  a  certain  value  of  a. 

The  results  discussed  above  for  Re(a40/aio)  and  Re(a4/a4)  differ  from  those  presented  in 
[5].  In  that  paper  it  was  reported  that  Re(a40/ai0)  and  Re(a4/a4)  were  always  negative,  so 
the  nonlinear  effects  were  always  stabilising.  Our  corrected  results  have  shown  the  significant 
result  that  the  nonlinear  effects  are  destabilising  for  particular  ranges  of  a.  Moreover  these 
ranges  correspond  to  realistic  values  of  the  location  of  porous  surfaces  in  experiments. 


8.2  Effect  of  the  shock 

Before  proceeding  to  investigate  the  effect  of  porous  walls  on  the  nonlinear  stability  it  is 
useful  to  consider  the  stability  problem  in  the  absence  of  shock.  The  problem  in  the  absence 
of  a  shock  was  first  considered  in  [12]  and  [14].  It  can  be  shown  that  the  neutral  curves  in 
the  absence  of  shock  differ  fundamentally  from  those  in  the  presence  of  shock  as  solutions  are 

34 


Distribution  A:  Approved  for  public  release;  distribution  is  unlimited. 


0 


-1.4  1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 - 1 

0.5  1  1.5  2  2.5  3  3.5  4  4.5  5 


(a) 


(b) 

Figure  23:  (a)  Re(a2o/al0);  (b)Re(a4o/aio)  for  the  first  four  modes  as  a  function  of  local 
cone  radius  a.  Results  are  shown  for  solid  wall  with  n  —  0  and  a/rs  =  0.57.  Results  for  the 
first  mode  are  indicated  by  o. 
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Figure  24:  Re(a4/ai)  as  a  function  of  local  cone  radius  a  in  the  absence  of  shock.  Results 
are  shown  for  solid  wall  with  (a)  n  —  1  and  (b)  n  =  2. 

only  possible  for  a  finite  range  of  a  (see  [4]).  If  no  shock  is  present,  solutions  to  the  upper- 
deck  equations  at  first  order  are  only  proportional  to  Kn(iar ),  allowing  only  for  outgoing 
waves  as  r  — >  oo.  This  modifies  the  resulting  eigenrelation.  The  nonlinear  stability  analysis 
can  be  carried  out  in  a  straightforward  manner  for  this  problem.  Figure  24  shows  Re(a4/ai) 
as  a  function  of  a  for  n  —  1  and  n  =  2.  It  can  be  seen  that  the  sign  of  Re(a4/ai)  is  always 
negative  and  two  solution  branches  exist  for  0  <  a  <  0.75  for  n  —  1  and  0  <  a  <  1.75  for 
n  —  2.  Thus  in  the  absence  of  shock  nonlinear  effects  are  stabilising  for  all  admissible  values 
of  a. 

8.3  Effect  of  porous  walls 

We  can  now  investigate  the  effect  of  porous  coatings.  We  begin  by  considering  the  regular 
porous  wall  model  of  (6).  The  results  using  the  regular  microstructure  model  comprising 
a  regular  array  of  cylindrical  pores  of  circular  cross-section  are  compared  to  the  results 
for  a  solid  wall  for  non-axisymmetric  modes.  The  flow  conditions  match  the  experimental 
conditions  of  [15].  The  porous  parameters  are  pore  radius  r*  =  28.5/im,  porosity  0O  =  0-2 
and  pore  depth  h*  =  450/irn.  Neutral  results  for  n  —  1  were  presented  in  figure  3a  which 
showed  that  the  porous  wall  reduces  the  neutral  values  of  The  porous  wall  has  only 
a  small  effect  on  the  neutral  values  of  a  (not  shown).  The  difference  between  the  neutral 
curves  of  the  solid  and  porous  walls  is  smaller  for  n  —  2  compared  to  n  —  1  especially  for 
the  lower  neutral  curves.  To  investigate  the  effect  of  nonlinearity  we  turn  to  figure  25  which 
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Figure  25:  Re(a4/ai)  for  the  first  five  modes  as  a  function  of  local  cone  radius  a  for  n  —  1 
and  a/rs  =  0.57:  — ,  solid  wall; - ,  regular  microstructure  model  (6). 

compares  Re(a4/ai)  for  n  —  1  between  the  solid  and  regular  porous  walls  for  the  first  five 
modes  (m  =  1,  5).  We  see  that  nonlinear  effects  are  enhanced  by  the  porous  wall  giving  larger 
values  of  Re(a4/ai).  In  the  presence  of  the  porous  wall  nonlinear  effects  destabilise  the  lower 
modes  (first  and  second)  while  stabilising  the  higher  modes  (three  to  five).  We  can  see  this 
from  the  fact  that  Re(a4/ai)  for  the  porous  wall  has  larger  positive  values  for  the  first  and 
second  modes  compared  to  the  solid  wall  and  that  Re(a4/ai)  for  the  porous  wall  becomes 
positive  at  larger  values  of  a  compared  with  the  solid  wall  for  the  higher  modes.  However  once 
destabilised,  the  higher  modes  of  the  porous  wall  have  larger  values  of  Re(a4/ai)  compared  to 
the  solid  wall.  For  large  enough  values  of  a,  we  can  expect  the  nonlinear  effects  to  diminish 
just  as  in  the  solid  wall  case.  Figure  26  shows  the  corresponding  results  for  Re(a4/ai) 
with  n  =  2.  Here  we  see  that  nonlinearity  destabilises  the  first,  second,  third  and  fourth 
modes  in  comparison  to  corresponding  modes  of  the  solid  wall,  while  stabilising  the  fifth  and 
possibly  higher  modes  of  the  porous  wall.  Results  obtained  for  n  —  3  (not  shown)  show  that 
nonlinearity  destabilizes  all  of  the  first  five  modes  of  the  porous  wall  compared  to  the  solid 
wall.  We  can  thus  infer  that  in  presence  of  the  porous  wall,  nonlinearity  destabilizes  lower 
modes  with  their  mode  number  increasing  with  azimuthal  wavenumber. 

We  now  compare  the  effect  of  porosity  on  the  nonlinear  stability  of  axisymmetric  distur¬ 
bances.  We  look  at  figure  27  which  shows  Re(a40/a1o)  as  a  function  of  a  for  the  porous  wall 
and  solid  wall.  Here  we  see  that  in  the  presence  of  porous  wall  all  the  modes  are  destabilised 
by  nonlinearity  with  the  most  significant  effect  being  felt  by  the  first  mode.  This  mode  is 
destabilised  for  all  values  of  a. 

Next  we  consider  the  random  microstructure  model  of  (11).  The  results  using  this  model 
are  compared  with  those  obtained  using  the  regular  microstructure  model.  The  porosity  of 
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Figure  26:  Re(a4/ai)  for  the  first  five  modes  as  a  function  of  local  cone  radius  a  for  n  —  2 
and  a/rs  =  0.57:  — ,  solid  wall; - ,  regular  microstructure  model  (6). 


Figure  27:  Re(a4o/aio)  for  the  first  four  modes  as  a  function  of  local  cone  radius  a  for  n  —  0 

and  a/rs  =  0.57:  — ,  solid  wall; - ,  regular  microstructure  model  (6).  Results  for  the 

first  mode  are  indicated  by  •  for  solid  wall  and  o  for  porous  wall. 
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Figure  28:  Re(a4/ai)  for  the  first  five  modes  as  a  function  of  local  cone  radius  a  for  n  —  1 

and  a/rs  =  0.57:  — ,  random  microstructure  model  (11); - ,  regular  microstructure 

model  (6). 

the  felt  metal  is  taken  to  be  0.75  and  the  fibre  diameter  is  30/mr.  For  comparison  the  regular 
microstructure  model  is  used  with  a  porosity  of  0.2  and  pore  radius  of  30/irn.  Figure  28 
shows  Re(a4/ai)  for  n  —  1  for  both  the  models.  We  see  that  nonlinear  effects  destabilise 
the  first  two  modes  of  the  felt  metal  model  when  compared  to  the  regular  porous  model. 
The  effect  on  the  higher  modes  is  opposite  as  we  see  that  Re(a4/ai)  becomes  positive  for 
smaller  values  of  a  for  the  regular  porous  model  as  compared  to  the  felt  metal  model.  Once 
destabilized,  values  of  Re(a4/ai)  are  more  positive  for  the  felt  metal  indicating  that  nonlinear 
amplification  of  disturbances  will  be  stronger.  Similar  trends  are  also  observed  (not  shown) 
for  the  case  of  axisymmetric  disturbances. 

Finally  we  consider  the  mesh  microstructure  model  of  (8).  The  results  using  this  model 
are  compared  with  those  obtained  using  the  regular  microstructure  model.  The  porosity 
of  the  mesh  model  is  0.8,  and  width  of  each  pore  section  is  100/irn.  For  comparison  the 
regular  porous  model  is  used  with  a  porosity  of  0.2  and  pore  radius  of  30/mr.  Figure  29 
shows  Re(a4/ai)  for  n  —  1  for  both  the  models.  The  effect  of  the  two  models  on  the  first 
mode  is  similar  with  the  destabilising  effects  of  nonlinearity  being  stronger  for  the  mesh 
model  compared  to  the  regular  model.  For  the  higher  modes  nonlinear  effects  are  slightly 
more  destabilising  for  the  regular  model  compared  to  the  mesh  model.  This  can  again  be 
seen  by  noting  that  Re(a4/ai)  becomes  positive  for  smaller  values  of  a  for  the  regular  model 
compared  to  the  mesh  model.  At  large  values  of  a,  Re(a4/ai)  is  more  positive  for  the  mesh 
model  for  the  first  four  modes  while  for  the  fifth  and  possibly  higher  modes  the  Re(a4/ai) 
is  larger  for  the  regular  model.  Corresponding  results  for  axisymmetric  disturbances  (not 
shown)  indicate  no  discernible  difference  between  the  effects  of  the  two  models. 
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Figure  29:  Re(a4/ai)  for  the  first  five  modes  as  a  function  of  local  cone  radius  a  for  n  —  1 

and  a/rs  =  0.57:  — ,  mesh  microstructure  model  (8); - ,  regular  microstructure  model 

(6). 

9  Discussion  and  conclusions 

We  have  presented  neutral  results  for  the  wavenumber  and  frequency  of  linear  disturbances 
for  hypersonic  flow  over  a  passive  porous  wall,  with  scales  appropriate  to  the  first  instability 
mode  for  axisymmetric  and  non-axisymmetric  disturbances.  The  porous  wall  parameters 
were  chosen  to  correspond  to  previous  theoretical  and  numerical  studies.  For  the  values 
chosen  here  we  find  that  the  neutral  values  of  streamwise  wavenumber  a  are  slightly  reduced 
from  those  corresponding  to  a  solid  wall.  However,  the  neutral  values  of  the  disturbance 
frequency  H  for  non-axisymmetric  disturbances  can  be  substantially  lower  than  those  for  a 
solid  wall,  particularly  for  larger  values  of  radius  a.  In  addition,  the  spatial  growth  rates 
presented  demonstrate  that  the  porous  wall  has  a  destabilising  effect  on  the  axisymmetric 
and  non-axisymmetric  modes.  The  growth  rates  are  significantly  larger  than  those  for  the 
solid  wall. 

The  effect  of  rarefaction  for  a  regular  microstructure  was  shown  to  be  destabilising. 
This  has  been  shown  by  [7]  to  give  larger  first  mode  destabilisation  and  less  second  mode 
stabilisation  for  hypersonic  flow  over  a  flat  plate. 

We  have  also  investigated  the  effect  of  the  porous  wall  parameters  on  the  maximum  spatial 
growth  rate.  The  effect  of  porosity,  pore  radius  and  pore  depth  on  the  maximum  spatial 
growth  rates  was  also  considered.  It  was  shown  that  increasing  each  of  these  quantities  leads 
to  larger  growth  rates,  but  they  level  off. 

We  have  made  some  comparisons  of  the  dimensional  growth  rates  for  the  first  mode 
disturbances  enhanced  by  the  porous  surface  with  those  obtained  from  previous  numerical 
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results  for  the  second  mode.  It  was  found  that  the  first  mode  growth  rates  were  smaller  than 
the  second  mode  ones. 

We  have  determined  the  effect  of  regular  microstructures,  random  microstructures  and 
mesh  microstructures  on  the  first  Mack  mode  in  a  hypersonic  boundary  layer.  Our  formula¬ 
tion  will  allow  for  alternative  porous  walls  to  be  investigated  in  a  straight-forward  manner. 

The  linear  stability  analysis  presented  will  not  be  valid  for  larger  disturbances.  Thus,  it  is 
important  to  determine  the  effect  of  nonlinearity  on  the  stability  of  hypersonic  boundary  layer 
flow  over  a  sharp  slender  cone  with  a  porous  wall.  This  has  been  investigated  experimentally 
for  second  mode  disturbances  using  the  bicoherence  method  in  [16]  and  [17]. 

The  weakly  nonlinear  stability  of  Mack’s  first-mode  disturbances  in  the  hypersonic  bound¬ 
ary  layer  on  a  sharp  slender  cone  with  passive  porous  walls  has  been  investigated.  The  effect 
of  the  attached  shock  is  found  to  be  significant.  In  the  absence  of  shock,  unstable  solutions 
are  possible  only  for  a  finite  range  of  cone  radius  a  and  nonlinearity  stabilises  linearly  unsta¬ 
ble  disturbances  for  all  admissible  values  of  a.  The  presence  of  the  shock  leads  to  multiple 
unstable  modes  for  all  values  of  a.  The  effect  of  nonlinearity  is  dependent  on  the  mode 
number  and  cone  radius  a.  For  axisymmetric  disturbances  on  a  solid  wall,  nonlinear  effects 
tend  to  stabilise  all  higher  modes  while  the  lowest  mode  is  slightly  destabilised  when  a  be¬ 
comes  large  enough.  In  the  presence  of  the  porous  wall  all  the  modes  are  destabilised  when 
compared  to  the  solid  wall.  Linear  stability  results  show  that  the  lowest  mode  is  the  most 
unstable  and  has  largest  spatial  growth  rates  for  both  solid  and  porous  walls.  This  most 
dangerous  mode  is  also  the  most  destabilized  by  nonlinearity  in  the  presence  of  the  porous 
wall.  For  non-axisymmetric  disturbances  on  a  solid  wall,  nonlinear  effects  destabilise  the 
lowest  mode,  while  the  higher  modes  are  stabilised  until  a  certain  value  of  a  which  increases 
with  azimuthal  wavenumber.  All  porous  wall  models  destabilise  the  neutral  modes.  When 
considering  the  effect  of  nonlinearity  on  linearly  unstable  modes,  we  can  state  that  lower 
modes  are  greatly  destabilised  by  nonlinearity  while  it  has  a  stabilising  effect  on  the  higher 
modes.  We  have  shown  that  it  is  the  higher  modes  that  have  the  largest  spatial  growth 
rates  in  the  presence  of  the  porous  wall.  The  effect  of  nonlinearity  is  to  stabilize  these  most 
linearly  amplified  modes  by  pushing  the  point  of  subcritical  instability  to  larger  values  of  a. 

The  random  microstructure  felt  metal  model  was  compared  with  the  regular  porous 
model.  The  felt  metal  significantly  destabilizes  the  neutral  modes  and  strongly  amplifies  the 
linearly  unstable  modes  with  the  higher  modes  giving  the  largest  growth  rates.  Nonlinear 
effects  in  the  presence  of  the  felt  metal  wall  stabilize  these  more  dangerous  higher  modes 
over  a  larger  range  of  a  while  destabilizing  the  more  slowly  growing  lower  modes.  When 
comparing  the  difference  between  the  mesh  microstructure  model  and  the  regular  porous 
model,  we  notice  similar  effects  between  both  models  on  all  the  modes.  The  regular  porous 
model  slightly  destabilises  all  modes  when  compared  to  the  mesh  model.  Since  the  felt  metal 
and  mesh  coating  have  higher  porosity,  to  corroborate  these  findings,  nonlinear  stability 
results  for  the  regular  porous  model  with  a  higher  porosity  of  do  =  f  were  obtained.  In 
figures  30a  and  30b  for  n  —  0  and  n  —  1,  respectively  we  see  that  higher  porosity  leads  to 
nonlinearity  having  a  stabilizing  effect  on  mode  numbers  greater  than  one.  This  is  a  result 
of  the  increase  in  the  value  of  a  where  Re(a4/ai)  becomes  positive.  However,  for  large  values 
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of  a  the  destabilising  effect  of  nonlinearity  is  stronger  with  Re(a4/ai)  being  slightly  larger 
for  higher  porosity.  In  figure  31a  for  n  =  2  we  see  the  stabilising  effect  for  mode  numbers 
greater  than  two  and  in  figure  31b  for  n  —  3  we  see  it  for  mode  numbers  greater  than  three. 
Thus  porous  coatings  with  higher  porosity  allows  nonlinear  effects  to  stabilise  higher  mode 
number  disturbances  at  a  particular  location  with  the  mode  number  of  the  lowest  mode  that 
is  stabilised  increasing  with  increasing  azimuthal  wavenumber. 

We  can  deduce  the  significance  of  curvature  on  the  nonlinear  stability  by  comparing  our 
results  for  a  solid  wall  to  those  obtained  from  the  analysis  of  [18]  for  the  weakly  nonlinear 
stability  of  flow  over  a  wedge.  Figure  32  shows  the  variation  of  the  Re(a3),  the  coefficient 
of  the  nonlinear  term  of  the  amplitude  equation  in  their  paper  (cf.  equation  5.1)  with  /3i, 
the  leading  order  scaled  spanwise  wavenumber.  The  results  are  shown  for  a  scaled  shock 
distance  ys  =  1.73  which  corresponds  to  the  shock  angle  expected  from  the  flow  conditions 
considered  in  this  paper.  Corresponding  results  (cf.  figure  5)  shown  in  [18]  were  obtained 
using  incorrect  values  for  two  constants  in  their  equation  4.12.  Our  corrected  results  shows 
the  significant  result  that  for  a  small  range  of  0  <  /3i  <  0.7,  Re(a3)  >  0.  Thus  nonlinear 
effects  will  be  destabilizing  for  disturbances  with  these  spanwise  wavenumbers.  By  comparing 
the  magnitudes  of  Re(a3)  and  Re(a4/oi)  we  can  infer  that  curvature  has  the  effect  of  making 
the  nonlinear  effects  stronger. 

We  now  discuss  the  implications  of  our  results  for  transition  to  turbulence.  The  weakly 
nonlinear  analysis  carried  out  here  shows  that  in  the  presence  of  porous  walls,  lower-frequency 
first  Mack  modes  are  destabilised  by  nonlinearity  while  higher- frequency  first  Mack  modes 
that  are  destabilised  on  the  solid  wall  at  a  particular  location  now  become  stabilised  for  a 
range  of  a.  This  effect  is  enhanced  by  models  with  higher  porosity.  The  porous  walls  amplify 
unstable  growth  rates  of  the  lower-frequency  modes  when  compared  to  the  solid  wall.  These 
faster  growing  modes  are  further  amplified  by  nonlinearity  and  may  contribute  to  the  early 
breakdown  of  laminar  flow.  Bicoherence  diagrams  from  the  experimental  investigation  of 

[22]  show  that  in  the  low-frequency  range  (fi,  f-2  <  100kHz)  nonlinear  processes  proceed 
more  intensely  on  the  porous  wall  compared  to  the  solid.  The  authors  suggest  that  this 
may  be  caused  by  the  growth  of  the  low-frequency  disturbance  amplitudes  due  to  surface 
roughness.  Bicoherence  measurements  of  [16]  have  also  identified  a  nonlinear  interaction 
that  is  associated  with  the  destabilised  first  Mack  mode.  The  authors  however  state  that  the 
nonlinear  interaction  is  too  weak  to  adversely  affect  the  performance  of  porous  walls.  These 
authors  only  considered  measurements  from  the  maximum  mass-flow  fluctuation  location 
while  in  [22]  measurements  were  taken  throughout  the  entire  thickness  of  the  boundary 
layer.  As  the  porous  wall  weakens  resonant  interactions  in  the  maximum  fluctuation  layer, 

[23]  show  that  nonlinear  interactions  above  and  below  the  layer  start  to  play  a  major  role. 
Nonlinear  interaction  between  vortex  (first-mode)  waves  and  filling  of  the  low-frequency 
vortex-mode  spectrum  in  the  presence  of  porous  walls  have  also  been  found  by  theoretical 
analysis  of  Gaponov  and  Terekhova  [24],  Simulations  by  De  Tullio  &  Sandham  [25]  of 
transition  over  a  flat  plate  in  the  presence  of  oblique  first  Mack  mode  show  that  the  first¬ 
mode  grows  faster  than  second  Mack  modes  and  drives  the  flow  directly  to  a  turbulent  state 
by  nonlinear  interactions  similar  to  the  oblique-transition  scenario  found  at  supersonic  Mach 
numbers  (see  [26]  and  [27]).  De  Tullio  and  Sandham  [25]  state  that  the  first  mode  regains 
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Figure  30:  Re(a4/ai)  as  a  function  of  local  cone  radius  a  for  a/rs  =  0.57  and  (a)  n  =  0;  (b) 

=  1.  Results  are  shown  using  the  regular  microstructure  model  (6):  — ,  0o  =  tt/4; - , 

)  =  0.2.  Symbols  refer  to  mode  number:  x,  m  =  1;  □,  m  =  2;  o,  m  =  3;  o,  m  =  4;  •, 
m  =  5. 
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Figure  31:  Re(a4/ai)  as  a  function  of  local  cone  radius  a  for  a/rs  =  0.57  and  (a)  n  =  2;  (b) 

=  3.  Results  are  shown  using  the  regular  microstructure  model  (6):  — ,  0o  =  7r/4; - , 

)  =  0.2.  Symbols  refer  to  mode  number:  x,  m  =  1;  □,  m  =  2;  o,  m  =  3;  o,  m  =  4;  •, 
m  =  5. 
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Pi 


Figure  32:  Re (0,3)  for  the  first  five  modes  as  a  function  of  spanwise  wavenumber  for 
ys  =  1.73  (cf.  equation  5.1  of  Seddougui  &  Bassom  1994). 

importance  in  the  transition  process  at  high  Mach  numbers  for  porous  surfaces.  The  results 
of  the  studies  reviewed  here  qualitatively  support  the  conclusions  of  our  analysis  that  the 
presence  of  porous  coatings  affects  the  instability  of  viscous  first-mode  disturbances  leading 
to  significant  nonlinear  development  of  disturbances  in  the  low-frequency  spectrum. 

From  this  discussion  it  is  clear  that  further  research  is  required  to  establish  whether 
Mack’s  first  mode  instability  can  cause  transition  to  turbulence  in  the  presence  of  porous 
walls.  The  work  presented  here  should  be  extended  to  consider  resonant  interactions  between 
the  unstable  modes  identified  in  this  study  and  its  harmonics. 

Since  it  has  been  observed  experimentally  and  confirmed  theoretically  and  numerically 
that  a  porous  microstructure  destabilises  the  first  Mack  mode,  it  would  be  desirable  to  design 
a  microstructure  which  stabilises  the  second  Mack  mode  while  not  appreciably  destabilising 
the  first  Mack  mode.  In  their  efforts  to  determine  types  of  microstructure  which  will  lead 
to  reduced  amplification  of  the  first  Mack  mode,  Wang  and  Zhong  [28]  investigated  the  ef¬ 
fect  of  the  phase  angle  of  the  wall  admittance  on  the  growth  of  the  first  Mack  mode.  For 
their  numerical  investigations  on  a  flat  plate  they  discovered  that  a  smaller  phase  angle  of 
the  wall  admittance  leads  to  a  weaker  destabilisation  of  the  first  Mack  mode.  The  effect 
of  wall  admittance  was  also  investigated  theoretically  by  Carpenter  and  Porter  [29].  They 
considered  a  thin  porous  sheet  stretched  over  a  plenum  chamber.  When  the  phase  of  the 
wall  admittance  was  very  close  to  n/2  they  discovered  that  Tollmien-Schlichting  waves  were 
completely  stabilised.  Thus,  future  studies  should  be  focused  on  determining  the  wall  pa¬ 
rameters  which  will  minimise  the  destabilising  effect  of  a  porous  wall  on  the  first  Mack  mode. 
This  will  enhance  the  current  knowledge  of  how  the  first  Mack  mode  in  hypersonic  boundary 
layers  may  be  controlled.  It  is  hoped  that  such  studies  will  provide  optimum  design  of  porous 
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coatings  for  incorporation  into  thermal  protection  systems  as  effective  techniques  for  laminar 
flow  control  in  hypersonic  flight  vehicles. 
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Appendix  A 

The  complex  constants  in  the  amplitude  equation  (16)  involve  the  following  expressions, 

where  as  defined  in  [19] 

K(0  =  Ai(£)  -  AlML(f),  (21a) 

where 

T(0  =  Ai(f)  f  /  Ai (q)dq,  (21b) 

Jto  Ai  (s)  J oo 

/**(£)  =  *1/3(Ai(0)(c)  {PiA^Ai(£o)  +  Ai'(0  -  Ai'teo)  -  £  f  M(t)dt\  +  c.c. 
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A12°  =  a-4/3T14A90A10°  +  a~4/3T 7  -  AyP10T13, 
A13°  =  T23  -  T24A9°A100, 

Pio  =  (Ay  +  ia)  1 , 

P20  =  (Ay  +  2ia)  , 

_  In(iars)Kn(iaa)  -  In(iaa)Kn(iars) 
In{ioirs)K'n{ioia)  —  I'n{iaa)Kn{iars)  ’ 

_  hn{2iars)K2n{2iar)  -  /2n(2iar)/i2n(2iars) 
5f2"  I'2n{2iaa)K2n{2iar s)  -  I2n{2iar s)K'2n(2iaa)  ’ 

0  _  /o(iar5)/i0(icrr)  -  /0(io;r)/io(mrs) 
n  I0(iars)K'0(iaa )  —  lQ(ioia)K'Q(ioirsy 
0  /o(2iars)A'o(2ia;r)  —  Jo(2iar)Ao(2ia;rs) 

^2n  I0(2iars)K'0(2iaa )  —  I'0(2iaa)  K'0(2iars) 
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Appendix  B 


The  resulting  amplitude  equation  for  axisymmetric  disturbances  for  fixed  frequency  is 

cL4.ii 


aio- 


dX 


—  020^2^11  +  fl40^4ll  |  ^4ll  |  • 


The  coefficients  in  (22)  may  be  expressed  as 
T1 

aio  =  -i4/3 - 2(ia)1/3Ai(^0) 

K 

(b0fo  +  c0d0)  (/lo(faa)  -  Al  —  (ia)~1/3F10K'0(iaa) 


,-i 

10 


x 


K 


( b0do  +  c0e0)  iloiiaa)  -  A*  ^  {ia)  1/3Floro(iaa) 


a2 o  =  i  2^olL‘1  +  an  lAYi1^3T3, 

r2/3  (a)5/3 


040  — 


K\K\ 


(22) 


(23) 

(24) 


ill"  -  Ay  [a2/3P10T22  -  P20X25LI20]  -  4v<c)Ptoil3°  ,  (25) 


where  the  corresponding  axisymmetric  versions  of  the  abbreviations  are  defined  in  Appendix 
A. 
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List  of  Symbols,  Abbreviations,  and  Acronyms 


a  = 

a  = 

&OC 

A 

Ay  — 

Ay 

Cd 

d* 

h 

L* 

Kn 

M 

n  = 

Pr  = 

P 
R 

Re  = 

rp 

rs  = 

T 

u,  v  = 

U 

x,r,4>  = 

a  = 

7 

ec 

6S 

A 

A 

to  = 

n 

V  = 

p 

Pd 

00  = 

Subscripts 


cone  radius 
half-pore  width 
tortuosity 

displacement  function 

admittance 

scaled  admittance 

dynamic  compressibility 
fibre  diameter 

h*/S *,  non-dimensional  porous-layer  thickness 

length  scale 

Knudsen  number 

Mach  number 

azimuthal  wavenumber 

Prandtl  number 

pressure  disturbance 

US* /v Reynolds  number 

Reynolds  number 

r*/8*,  non-dimensional  pore  radius 

shock  location 

temperature 

velocity  disturbance 

velocity 

orthogonal  coordinates 

characteristic  impedance 

ar  +  disturbance  streamwise  wavenumber 

boundary-layer  displacement  thickness 

specific  heat  ratio 

cone  angle 

shock  angle 

skin  friction 

propagation  constant 

angular  frequency 

frequency 

viscosity 

kinematic  viscosity 
density 

P*d/Pw i  dynamic  density 

porosity 

flow  resistivity 
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—  =  just  behind  the  shock 

s  =  shock 
w  =  wall 
0  =  initial 
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